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^ \ Abstract: We estimate the current theoretical uncertainty in sparticle mass pre- 

Q-i! dictions by comparing several state-of-the-art computations within the minimal su- 

Qh! persymmetric standard model (MSSM). We find that the theoretical uncertainty is 

r-| ! comparable to the expected statistical errors from the Large Hadron Collider (LHC), 

and significantly larger than those expected from a future e^e~ Linear Collider (LC). 
^ \ We quantify the theoretical uncertainty on relevant sparticle observables for both 

^ \ LHC and LC, and show that the value of the error is significantly dependent upon 

the supersymmetry (SUSY) breaking parameters. We also present the theoretical 
uncertainty induced in fundamental scale SUSY breaking parameters when they are 
fitted from LHC measurements. Two regions of the SUSY parameter space where 
accurate predictions are particularly difficult are examined in detail: the large tan j3 
and focus point regimes. 
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1. Introduction 

Weak-scale supersymmetry (SUSY) is well motivated because it solves the 
technical hierarchy problem, removing ultra-violet quadratic divergent corrections 
to the Higgs mass. Viable weak-scale supersymmetry implies over a hundred "soft 
breaking" terms in the MSSM Lagrangian at the weak scale, making general analysis 
intractable. By making model-dependent assumptions on the origin of SUSY break- 
ing for the MSSM fields, one often provides relations between the weak-scale SUSY 
breaking parameters, vastly reducing the available parameter space. Many theoreti- 
cal schemes of SUSY breaking exist, the most popular of which can be classified by 
the mechanism that mediates SUSY breaking from a hidden sector to the MSSM 
fields. In this paper, we will use minimal supergravity (mSUGRA) minimal 
anomaly mediation (mAMSB) ^ and minimal gauge (mGMSB) [^] mediation. 

Supposing SUSY is discovered at a present or future collider, it will be a major 
challenge to measure the SUSY breaking parameters with good accuracy. In order 
to determine the free parameters, physical sparticle masses (or kinematical variables 
related to them) will be measured. These sparticle masses must be turned into SUSY 
breaking parameters, typically at some high scale. Provided enough information is 
collected, this will allow tests on the relations between the high-scale parameters. For 
example, a simple question already addressed [§, 0; is: do the gaugino masses 
unify, and if so, do they unify at the same scale as the gauge couplings? A double 
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affirmative would be strong support in favour of SUSY GUT-type schemes. On the 
other hand, there are many other possibihties: for example, string models can be 
non-universal |]10|, or unify at an intermediate scale. 

There is a large literature upon expected empirical errors on the observables 
both a Linear Collider 1T1|, |12[ |T3|, [T|, g and the LHC |16|, [l^, 0, ig. It will 



be essential to know the theoretical uncertainty if we are to discriminate models of 
SUSY breaking. There are several studies of how empirical errors propagate into 



uncertainties on fundamental-scale SUSY breaking parameters in the literature pO 
|l],||,|l6|,0,|3|,gi. 



The deduction of high-scale SUSY breaking parameters from observables in- 
evitably involves theoretical errors coming from the level of approximation used (e.g. 
neglected higher order terms), and it is these uncertainties that we study here. We use 



four modern codes to calculate MSSM spectra: ISAJET7.64 gj, S0FTSUSY1.71 [gS 
SPHENO2.0 ra and SUSPECT2.101 [Bl. We use the differences in results between 



the codes to define the current theoretical uncertainties. We have done our best to 
eliminate differences due to bugs by examining the relevant parts of codes in detail 
if there was an obvious large discrepancy. However, it would be unrealistic to claim 
that all of the codes are completely bug free. We therefore take a practical inter- 
pretation of 'theoretical uncertainty': after all, when ffiting experimental data to 
SUSY breaking models, one must use one of the available computational tools. The 
precise implementation of the known higher-order corrections differs and has been 
found to produce significantly different results (for example, using different scales for 
parameters in the highest-order corrections). Therefore, certainly the differences in 
results between the codes are due (at least to a large part) to unknown higher-order 
corrections, matching more traditional notions of 'theoretical uncertainty'. 

Previously, there have been some studies of differences between various calcula- 



tions of the MSSM spectrum. 2-5% differences in sparticle masses were noticed ^8 



between various particles along the Snowmass (SPS) p9| model lines la and 6 be- 
tween the SOFTSUSY 1.2, ISAJET 7.51 and SUSPECT 2.0 programs. Model line 2 (the 
so-called focus point line, with very heavy scalar sparticles) was observed to show 
huge 30% differences in the masses of the weak gauginos. Differences in the spec- 
tra and branching ratios were also observed in comparing ISAJET 7.58, SUSY- 
GEN 3.00 and PYTHIA6.2. Moreover, large 10% level differences between SOFT- 
SUSY 1.3, ISAJET 7.58 and feynSSG were observed |T§ in the mSUGRA Post-LEP 
benchmarks [31|. Some of the benchmarks in some of the codes were found to not be 
consistent with electroweak symmetry breaking, unless one fiddled with the Standard 
Model inputs and as{Mz). These were noticeably points E,F (focus-points) and 
K,M (high tan/5 and high mo,mi/2 points). Likewise, the input parameters had to 
be adjusted in to get similar spectra from ISAJET 7.51 and SSARD. Initial results 
highlighting the differences between the predictions in the focus-point and high tan (3 
regimes have already been presented by the current authors as a conference proceed- 
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ing . We will present the main results of this last work here for completeness, 
updated to state-of-the art calculations. 

In this paper, we push these initial observations of theoretical uncertainties fur- 
ther by (i) comparing them to the expected experimental accuracies at the LHC and a 
future LC and (ii) examining their effect upon future empirical fits to fundamental- 
scale SUSY breaking schemes. We take expected statistical errors upon sparticle 
observables from previous studies of mSUGRA at the LHC [|T^] at two of the LHC 
benchmark points. We then perform a fit to mSUGRA at each point with each of the 
four codes. The statistical precision of the fitted fundamental mSUGRA parameters 
may then be compared with the theoretical error by looking at the differences be- 
tween the four fits. We next quantify and present the theoretical error in the coloured 
sparticle masses at the SPS points. We also re-examine the SPS points providing 
theoretical uncertainties on mass predictions for sparticles that may kinematically 
be accessed at a 500 GeV LC. Quantification of errors at particular points will give 
us an idea of their magnitude and whether or not they significantly depend upon the 
SUSY breaking scheme (or point). 

In section |^, we introduce the four state-of-the art calculations we will use and 
their level of approximation: ISAJET7.64, S0FTSUSY1.71, SPHENO2.0 and SUS- 
PECT 2.101. 'Tricky' regions of the MSSM parameter space, where it is difficult to 
make accurate predictions, are discussed in section ^ large tan P and the focus- 
point regime. In section [4.1| , we perform the LHC mSUGRA empirical fits and then 
quantify the theoretical error on squark and gluino masses for the SPS points. In 
section |4.2|, we quantify theoretical errors on masses relevant for a 500 GeV Linear 



Collider. Finally, there are conclusions and an outlook in section ^. 



2. The codes 

We compare the latest versions of four public SUSY renormalisation group evolution 
(RGE) codes which we think constitute a representative sample of such programs: 
ISAJET7.64, S0FTSUSY1.71, SPHENO2.0 and SUSPECT 2.101. The basic principle 
of the SUSY mass spectrum calculation is the same in all programs: Gauge and 
Yukawa couplings are taken as input parameters at the electroweak scale. However, in 
the MSSM, in order to define them in the DR scheme, from experimental observables, 
one must first subtract threshold corrections from sparticles. The sparticle spectrum 
is unknown at this stage and so to begin the calculation, some guess is made for the 
soft SUSY breaking parameters and spectrum. The MSSM parameters are then run 
to the high scale Mx by RGEs. At Mx, boundary conditions are imposed on the 
SUSY breaking parameters. Couplings and SUSY parameters are then run back down 
to MsusY = y/^Wj^ where the and B MSSM Higgs potential parameters are set 
in order to give correct radiative electroweak symmetry breaking consistent with an 
input value of tan/5 = V2/V1 (fi, V2 being the two Higgs fields' vacuum expectation 
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full 1-loop 


full 1-loop 
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full 1-loop 


full 1-loop 
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all have 1-loop gg + 


qq re-summed 





Table 1: RGEs and radiative corrections implemented in ISAJET 7.64, SUSPECT 2.101, 
S0FTSUSY1.71 and SPHENO2.0. 



values (VEVs)). The SUSY masses are calculated and radiative corrections are 
applied, and the parameters are run down to the electroweak scale. Finally, the 
whole process is iterated in order to obtain a stable solution. 

An overview of which corrections are implemented in each of the four programs 
is given in Table |1|. From this table we already expect some differences in the results 
due to the different levels of radiative corrections applied. In particular ISAJET 7.64 
has no finite radiative corrections to most of the sparticle masses. However, we 
note that even if each column of the table were identical, one could expect different 
numerical results from each of the codes. In practise, if a quantity is calculated 
at one-loop, one has the freedom (at one-loop accuracy) of using whatever scale one 
desires for parameters in the one-loop correction itself. The difference between using, 
for instance, pole or running masses or couplings derived in the MS or DR scheme 
in one-loop corrections is formally of higher-loop order, but leads to non-negligibly 
different results. Indeed, one can roughly estimate the effects of higher-loop terms by 
the difference in predicted masses by varying the scale of parameters in the highest 
loop included. 

Let us now point out some of the differences in more detail. A subtle but 
important issue is the treatment of Yukawa couplings. In general, they are derived 
from the quark masses as 

ht = V2mt/v2, hb = V2mb/vi, (2.1) 

where rrit and rrih are the running top and bottom quark masses in the DR scheme. 
Obviously, it makes a difference whether the VEVs vi^2 are running (S0FTSUSY1.71, 
SPHENO2.0, SUSPECT2.101) or not (ISAJET 7.64). It makes an important difference 
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how the DR masses are calculated from the pole or MS masses. We first discuss 
Standard Model threshold corrections, and afterwards the sparticle loop corrections. 
• ISAJET7.64 takes 2-loop QCD corrections to rrit into account, including the shift 



from the MS to the DR scheme as Il39 

5a. 



mtiM, 



DR 
t)SM 



1 + 



3n 



16.11 - 1.04 5 



6.63 



(2.2) 



where Mf is the top pole mass and = ^^(Mj) at 3-loops in the MS scheme. 
For the bottom quark mass, ISAJET7.64 takes a hard-coded value of mb{Mz)sM = 



2.82 GeV H. 

• S0FTSUSY1.71 and SPHENO2.0 calculate both ht and h at Q = Mz- The DR 
top mass is related to the pole mass by 2-loop QCD p 



miQ) 



DR 

SM 



3L) -af 0.538 - 



43 

247r2 



87r2 



(2.3) 



where^ L = \n{M^ /Q"^). For the h quark, 3-loop relations |^ are used to calcu- 
late mft(Mb)g^ from the h pole mass, which is run to Mz by 3-loop RGBs ^]. 
mi,{Mz)^^ is then related to the DR mass by |2^ 



MMz)fl 



37r 



+ 



727r2 1287r2 11527r2 



(2.4) 



in S0FTSUSY1.71. SPHENO2.0 neglects the last two terms of eq. (|J). 

• SUSPECT 2.101 uses 2-loop relations and 2-loop RGBs to derive mi,{Mb)f^ and 

mt(Mt)g5^ from the h and t pole masses. For the conversion to the DR scheme, 

eqs. ( p^.3| ) and ( p^.4D are applied (with L = in eq. ( P73| ) since Q = Mt). 

Using Q = Mt to define mt{Q) (ISAJET7.64, SUSPECT 2.101) is in principle more 

accurate than using Q = Mz (S0FTSUSY1.71, SPHENO2.0) since then the QCD 

logs between Mz and Mt are re-summed. We refer the reader to the respective 



manuals |2J, |2^, |2g, ^ for more details. 

The next step is to include SUSY loop corrections. For m^, all four programs ap- 



ply full 1-loop SUSY corrections according to |Q. For m^, the full 1-loop (ISAJET 7.64, 
SPHENO2.0) or the leading (S0FTSUSY1.71, SUSPECT 2.101) SUSY corrections are 
included. The tan/3 enhanced corrections to are re-summed as given in in 



^After the publication of this paper, it was noted that eq. (2_^) incorrectly describes the 2-loop 
QCD corrections to the top mass: eq. (2^) holds for L = \n{m.^{Q)/Q^), not for L = ln{Mi /Q^). 
S0FTSUSY1.71 and SPHENO2.0 therefore contain the incorrect formula, which ought to be 



3tt 



{5-SL)-ai[ 0.876 



91 
247r2 



87r2 



■L' 



(2.3a) 



for L ~ \n{M^ /Q'^). We thank D. R. T. Jones for drawing our attention to this issue. 
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all four programs, see [Q. Still, there are some important differences in the imple- 
mentation of these corrections. For example, the gt correction to rrit is 



(Arrit \^ { -rj / 2 2 2\iD/2 2 2\ 



+ ^ sin 29^ [Bo{ml mf, m|) - 5o(m?, m^, m|)] |. (2.5) 

IS A JET 7.64 adds this correction at m^, but using a scale Q = ^m^^mi^, (this is 
done to avoid double counting of logarithmic corrections which are included via step 
functions in the RGBs) and = as{mg) in eq. (PT^), while S0FTSUSY1.71 and 
SPHENO2.0 calculate it at Q = Mz and SUSPECT 2.101 at Q = Mj. Accordingly, 
nit = '^i(Q)sM the term rrig/mt enters with different values in all four programs. 
Due to differences in the inclusion of finite radiative corrections to sparticle masses, 
the gluino and stop masses in (|2.5|) vary from program to program. In particular, 
ISAJET7.64 calculates the corrections to nig aX Q = m-g while SUSPECT 2.101 cal- 
culates them at Q = M^, which leads to quite different gluino masses. Analogous 
differences exist in the other contributions to (Amf)"^^"^^ as well as in the calculation 
of (Amb)^^^^. 

Another comment is in order concerning a^. The value of as{Mz) from experi- 
ment is given in the MS scheme. S0FTSUSY1.71, SPHENO2.0 and SUSPECT 2.101 



take as{Mz) as input and convert it to the DR scheme |^ 



,r-iMz).'^^, (2.6, 



where 

Ois{Mz) 



27r 



Aa, 



1 2 , / Mf \ , / nia \ 1 \ ^ X ^ , / nia 



In — - 2 In 



2 3 \Mz) \MzJ Q^^, \Mz 

q 1=1,2 



-Win 



:2.7) 



In SUSPECT 2.101, the log terms are not added explicitly but included via threshold 
functions in the RGBs. Also ISAJET 7.64 re-sums the logs in ( |2.7| ) via step-function 
decoupling in the RGBs. The finite term, however, is not taken into account. The 
difference due to the finite term is small but relevant (1%). One could in princi- 
ple interpret the input as{Mz) in ISAJET 7.64 as already being the effective Stan- 
dard Model DR value. Since, however, for some corrections ISAJET 7.64 effectively 
takes hard-wired values of as{Mz)^^ = 0.118, we take the canonical input value of 
a^iMz)^ = 0.118 for all codes. 

3. Tricky corners of SUSY parameter space 

As a general point, when quantifying errors on predicted observables, we assume 
that the results from the codes follow a Gaussian probability distribution. This is a 
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0. 



10 20 30 40 50 

tan/? 

Figure 1: Bottom Yukawa coupling at MsusY scale as a function of tan /3, for tuq = 
400 GeV, mi/2 = 300 GeV, = 0, ^ > 0, Mt = 175 GeV; full (dotted) lines: 
ISAJET7.64 (7.58), dashed: SOFTSUSY 1.71, dash-dotted: SPHENO2.0, dash-dot- 
dotted: SUSPECT 2.101. 



priori unjustified, but we prefer it to quoting minimum and maximum values because 
we believe that the additional information included in the variance is desirable. For 
example, if three codes all provided identical results and one gave a result further 
away, we think that the true uncertainty ought to be less than the range of minimum- 
maximum results. 



3.1 Large tan (3 

Large tan (3 has always been recognised as a difficult case since it requires a thorough 
treatment of the bottom Yukawa coupling. Figure shows of the four different pro- 
grams as a function of tan (3 in the mSUGRA model. We see that SOFTSUSY 1.71 and 
SPHENO 2.0 agree very well on hi,, and there is also good agreement with ISAJET 7.64. 
Comparing only these programs we would assign a < 3% uncertainty on hi, even 
for very large tan/5. The agreement with SUSPECT 2.101 is however not so good, 
and we find 4-8% uncertainty taking all four programs into account. The effect of 
re-summing the tan /3 enhanced SUSY loop corrections can be seen when comparing 
the solid and dotted lines in fig. |l], the solid line being the result of ISAJET 7.64, 
where the (Am;,)'^^'^^ corrections are resummed, and the dotted one being the result 
of ISASUSY7.58, where this re-summation is not applied. 

The bottom Yukawa coupling has its largest effect in the Higgs sector when it is 
large (at high tan/3): the evolution of m|^^ is driven by /i^, 

dm? 3 

-^^—h,X, + ..., X, = {ml + ml + ml^ + Al), (3.1) 
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Figure 2: Higgs boson masses as a function of tan for rriQ = 400 Gev, mii2 = 300 GeV, 
Aq = Q, ^l > Mt = 175 GeV; full (dotted) lines: ISAJET7.64 (7.58), dashed: SOFT- 
SUSY 1.71, dash-dotted: SPHENO2.0, dash-dot-dotted: SUSPECT 2.101 (for the grey 
dash-dot-dotted hne corresponds to SUSPECT 2.101 + FeynHiggsFast) . 



where t = InQ, Q being the renormalisation scale. Differences in directly 
translate into the physical Higgs boson masses since 

= - - -?.J + ^ + ^ - Ml . (3.2) 

C2/3 Vi V2 

Here m|^. = mjj. — ti/vi, i = 1,2, and ti^2 are the tadpole contributions. The self 
energies of Z and A have been neglected in eq. (|3.2|) . 

Figure |^ shows the Higgs boson masses obtained by the four programs as a 
function of tan/5 (the results obtained by ISAJET7.58 are again shown as dotted 
lines in fig. 0). Let us first discuss the masses of A° and H^. For tan/? = 10-50, we 
find differences in and mfj± of about 10-50 GeV, dominated over most of the 
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tan /3 range by the difference between SPHENO2.0 and SUSPECT 2.101. This has to 
be compared with differences of 100 GeV and more encountered with earher versions, 
see for instance the dotted hnes representing ISAJET7.58. Assuming the error to be 
Gaussian, we now have /S.mA,H± — ±10 GeV at tan /5 = 25 and /S.mA,H± — ±20 GeV 
at tan/? = 50. The bottom Yukawa couphng is, however, not the only source of 
differences in m^. Another source is, for example, whether one uses running or pole 
values for masses in the calculation of one-loop tadpoles ti 2- Also, the scale and 
scheme of parameters in the one-loop expressions for the tadpoles all vary. These 
differences are formally of higher order and indeed each program has a different 
approach. 

The situation is somewhat different for the H^) system because here addi- 
tional radiative corrections are necessary. It is well known that these involve a theo- 
retical uncertainty on m/jo of about 3 GeV |^6|, evidence of which can be seen in fig. ^ 



For completeness we note that SUSPECT 2.101 offers various choices of Higgs mass 
calculations. In fig. |], we have used its default m^o routine, i.e. ichoice(lO) =0, 
shown as a black dash-dot-dotted line. If we use instead SUSPECT 2.101 with Feyn- 
HiggsFast, ichoice(lO) = 3, we get m/jO ~ 115 GeV, shown as a grey dash-dot-dotted 
line in fig. 0. This will be relevant later in this paper when we discuss mSUGRA fits 
to LHC data. 

3.2 Focus point 

For large mg, the running of m^^^ becomes very steep and very sensitive to the top 
Yukawa coupling: 

-^^—h,X, + ..., X, = {ml + ml + ml^+Al). (3.3) 
As a result, the parameter given by 

becomes extremely sensitive to ht. This is visualised in fig. ^ where we show ht and 
/i as functions of itiq. The other parameters are mi/2 = 300 GeV, = 0, tan/? = 10 
and /i > 0. There is reasonable agreement on /i up to mo ~ 1 TeV, and the differences 
observed for mo ^ 2 TeV are phenomenologically not so important. For larger values 
of mo we observe, however, large discrepancies between the four programs. These 
lead to completely different chargino/neutralino properties for very large mo, and 
likewise to very different excluded regions, depending on which program is used. 
The situation is, however, already much better than the one reported in c.f. the 
dashed lines of ISASUSY7.58 in fig. ^ (earlier versions of the other codes also gave 
quite different results). 
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Figure 3: ht and /x as a function of mo for = 300 GeV, ^o = 0, tan/? = 10, 

/i > 0, Mt = 175 GeV; full (dotted) lines: ISAJET7.64 (7.58), dashed: S0FTSUSY1.71, 
dash-dotted: SPHENO2.0, dash-dot-dotted: SUSPECT 2.101. 
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Figure 4: (a) as given by eq. (|3.5[ ) for mi/2 = 300 GeV, = 0, and tan/3 = 10 in 
the {mo,rht) plane, (b) Gluino mass as a function of mo for mi/2 = 300 GeV, Aq = 0, 
tan/3 = 10, ^ > 0, Mt = 175 GeV; full lines: ISAJET7.64 dashed: S0FTSUSY1.71, 
dash-dotted: SPHENO2.0, dash-dot-dotted: SUSPECT 2.101. 



In order to understand the behaviour of /i in fig. |] it is useful to write eq. (|3] 
in the form 

/x^ ~ cimo + C2mi/2 -0.5M|. (3.5) 
Approximate analytical expressions for ci and C2 can be found e.g., in [E^, E3]. For 
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Ao = and tan /5 = 10 we get [|| 



!] 



V 156.5 GeV/ ' V 102.5 GeV/ ^ ^ 

Figure ^ shows contours of constant fi in this approximation in the (mo, rrit) plane. 
Notice the fast increasing dependence on for increasing mo. For m^ ~ 156- 
157 GeV, fi becomes almost independent of mo. This is a signal of the actual focus 
point behaviour [^, which is defined as the insensitivity of m^/j to its GUT scale 
value. Eq. |3.4| then implies that fi is insensitive to mo, since appears with 
a suppression factor of tan^/5 — 1. Interpreting m^ in eq. ( p.6|) as mt{MsusY) = 
ht{MsusY) V2{Msusy) / we can directly relate the mo dependence of n in fig. ^ 
to that of ht- 

Some more comments are in order. Firstly, SUSPECT 2.101 has only 1-loop 
RGEs for scalar SUSY parameters. For large mo, the 2-loop terms lead to (9(10%) 
correction and should thus be taken into account. Secondly, as already mentioned 
in section 0, the sparticle masses that enter the radiative corrections have different 
values in different codes. In particular, there are large differences in the gluino 
masses between SUSPECT 2.101 and the other codes, as illustrated in fig. |b. Since 
the gluino mass enters ht via the gt correction eq. (|2.5| ), this may account for the 
different slope of hf (and consequently of /i) in SUSPECT 2.101 compared to the other 
programs, as evident in fig. 0. 

It thus seems that for reliable results in the large mg region, a more complete 
calculation of the top Yukawa coupling at the 2-loop level is necessary. This should 
include finite radiative corrections to all sparticle masses at the full 1-loop level. 



4. Comparison of theoretical and experimental uncertainties 
4.1 Fits of mSUGRA parameters to LHC data 

In the ATLAS TDR [|l|], study was made of fitting the mSUGRA model 



to possible measurements of six reference scenarios. We have re-analysed these fits 
for two of these scenarios, LHC Point 1 and Point 2. Here squarks and gluinos are 
produced with the dominant decays g qq^^ j^, ^ Xi^q, q^ Xil- The 

assumed measurements for the two points for low {C = 30fb~^) and high {C = 
300 fb~^) luminosity are given in table ^ They were estimated in ref. [|1^ by using 



ISAJET 7.34 and simulating the ATLAS experiment to determine expected empirical 
errors.^ 

With each of the programs under discussion we have performed a fit of 
mSUGRA parameters mo, mi/2 and tan/5 to the data of table |^, taking Aq = 

^It is beyond the scope of the present paper to re-perform the experimental analysis in order to 
have the numbers in table more up-to-date. 
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and n > 0. The results are listed in tables and ^. The quoted errors are at 
la (68.3% CL.) from a simultaneous fit of all three parameters, i.e. Ax^ = 3.53. 
In case of SUSPECT2.101, we have used its default option, ichoice(lO) =0, for 
the calculation of the mass. When linking SUSPECT 2.101 with FeynHiggsFast, 
ichoice(lO) =3, the results for mo and mi/2 practically do not change. However, we 
get much lower values for tan/3: tan/3 = 1.7±0.1 for Point 1, and tan/3 = 6.06±2.06 
for Point 2 at high luminosity. 

It is interesting to note that not only the central values but also the size of the 
errors can be quite different. Similarly, the minimum x^, which is a measure of the 
quality of the fit, can show large variations. To make the comparison easier and to 
visualise correlations and non-Gaussian effects, we show in Figs. || and || contours of 
68% and 95% CL. for the fits of Points 1 and 2 in the (mo, mi/2) and (mo, tan/3) 
planes. The values of Ax^ used for these confidence levels are based upon a simul- 
taneous two-parameter fit, i.e. Ax^ = 2.3 and 5.99. The third parameter, tan/3 or 
mi/2, is always fixed to its best-fit value. As one can see, the error ellipses have only 
little or even no overlap. We therefore conclude that the theoretical uncertainty is 
about the same size as the statistical one in the fitted quantities. 

One might expect that the main source of these differences is the theoretical 
uncertainty on m/j. This is indeed the case for the determination of tan/3. However, 
rrih has only little influence on the flt of mo and mi/2. This becomes clear when 
using e.g., SUSPECT 2.101 with different routines for the Higgs mass calculation. In 
flg. |a, we show the results of SUSPECT 2.101 FeynHiggsFast, ichoice(lO) =3, as 
dashed contours in addition to those obtained with its default Higgs mass routine 
(ichoice(lO) = 0, solid contours). In flg. Pd, we have omitted the default SUS- 
PECT 2.101 results to avoid confusion of the many lines. They would look similar 
to the SPHENO2.0 contours but centred at mo = 450 GeV and without an upper 
limit on tan/3. We thus conclude that the uncertainties in mo and mi/2 mainly come 
from the differences in the programs pointed out in Sects. 0-^, and not from different 
Higgs mass calculations. 

We next address the question of how the theoretical uncertainty depends on the 



Quantity 


Low-L 


High-L 


ruh (Point 1) 


95.4 ± 1.0 GeV 


95.4 ± 1.0 GeV 


rrih (Point 2) 


115.3 ±1.0 GeV 


115.3 ± 1.0 GeV 


^hq 


758.3 ± 25 GeV 


758.3 ± 25 GeV 




959 ± 40 GeV 


959 ± 15 GeV 




1004 ± 25 GeV 


1004 ± 12 GeV 


(Point 1) 


none 


647 ± 100 GeV 


mi^ (Point 2) 


none 


713 ± 100 GeV 



Table 2: Possible LHC measurements for Point 1 and Point 2, from [16| 



Point 1, Low-L 





• 

' ^ 1 1 LLt L 


mo [GeV] 


mi/2 [GeV] 


tan (3 


ISAJET 


0.10 


490 ± 135 


424 ± 25 


1.77 ±0.21 


SOFTSUSY 


9.30 


280 ± 246 


425 ± 32 


1.60 ±0.03 


SPHENO 


0.02 


373 ± 175 


436 ± 26 


2.10 ±0.15 


SUSPECT 


0.32 


411 ± 116 


410 ±20 


2.08 ±0.16 






Point 1, Hij 








• 


mo [GeV] 


nil 12 [GeV] 


tan (3 


ISAJET 


0.57 


496 ± 61 


424 ± 12 


1.77 ±0.20 


SOFTSUSY 


11.66 


356 ± 78 


422 ± 12 


1.60 ±0.03 


SPHENO 


0.27 


370 ± 82 


436 ± 12 


2.10 ±0.15 


SUSPECT 


1.79 


422 ± 67 


409 ± 13 


2.08 ±0.15 


Table 3: Fit to LHC measurements for Point 1. 






Point 2, Low-L 






Amm 


mo [GeV] 


mi/2 [GeV] 


tan /? 


ISAJET 


0.03 


523 ± 129 


424 ± 24 


6.55 ±2.37 


SOFTSUSY 


0.08 


414 ± 140 


419 ±23 


4.65 ±0.76 


SPHENO 


0.19 


405 ± 167 


437 ± 26 


> 7 


SUSPECT 


0.06 


444 ± 114 


409 ± 18 


> 7 






Point 2, Hij 


5h-L 






A.min 


mo [GeV] 


mi/2 [GeV] 


tan (3 


ISAJET 


0.14 


521 ± 58 


424 ± 11 


6.52 ±2.30 


SOFTSUSY 


0.33 


411 ±68 


419 ± 11 


4.63 ±0.88 


SPHENO 


1.08 


394 ± 80 


438 ± 13 


> 7 


SUSPECT 


0.20 


450 ± 64 


408 ± 11 


> 7 



Table 4: Fit to LHC measurements for Point 2. 



SUSY parameter point. Copious quantities of squarks and gluinos are expected to 
be produced at the LHC, leading to a fairly precise measurement of their masses, 
particularly if m^o is determined accurately by a LC. In table ^ we compare the g, m^. 



Uj^ and ti masses obtained by the four programs for the Snowmass (SPS) points [^ 



Assuming a Gaussian distribution, we quote the variance of these masses as the 



theoretical error, i.e. 5mx = y j^'^i[{^x)i — rnx]"^ where mx is the mean of 
{mx)i and = 4 in our case. We make the following observations: (i) the absolute 
theoretical uncertainty (in GeV) varies from point to point; (ii) the typical relative 
uncertainty in mSUGRA and mGMSB scenarios in generic (i.e. not tricky) regions 
of parameter space is about 2-5%; (iii) in some cases, in particular in focus point 
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Figure 5: Fit to LHC measurements for Point 1 with ISAJET 7.64, S0FTSUSY1.71, 
SPHENO 2.0 and SUSPECT 2.101 for £ = 300 fb'^ Shown are contours of 68% and 95% 
C.L. for each program in the {mo, mi/2) and (mo, tan/3) planes, with the third parameter 
fixed to its best fit value, c.f. table U. 
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Figure 6: Fit to LHC measurements for Point 2 with ISAJET 7.64, SOFTSUSY 1.71, 
SPHENO 2.0 and SUSPECT 2.101 for C = 300 fb'^ Shown are contours of 68% and 95% 
C.L. for each program in the {mo, mx/2) and ("lO) tan/3) planes, with the third parameter 
fixed to its best fit value, c.f. table In case of SUSPECT 2.101, the solid lines are for its 
default m/i routine and the dashed lines for mh calculated with FeynHiggsFast. 



and mAMSB scenarios, the relative uncertainty is larger, about 5-10%; (iv) in any 
case, the theoretical error is of the same order of magnitude as the experimental one. 



4.2 Linear Collider measurements 

At a high-luminosity e^e~ Linear Collider, one expects to measure chargino, neu- 
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mass 


code 


la 


lb 


2 


3 


4 


5 


6 


7 


8 


9 


g 


ISAJET 


607 


936 


794 


932 


732 


719 


718 


944 


835 


1296 




SOFTSUSY 


614 


949 


802 


946 


743 


730 


729 


964 


852 


1306 




SPHENO 


594 


917 


782 


914 


719 


705 


704 


940 


836 


1232 




SUSPECT 


626 


964 


870 


959 


761 


730 


742 


986 


902 


1395 




error 


13 


20 


40 


19 


18 


12 


16 


22 


31 


67 


Ur 

Li 


ISAJET 


536 


835 


1532 


817 


730 


642 


640 


858 


1079 


1233 




SOFTSUSY 


549 


851 


1582 


831 


753 


657 


662 


876 


1083 


1291 




SPHENO 


565 


876 


1563 


859 


764 


676 


674 


910 


1127 


1314 




SUSPECT 


570 


886 


1595 


867 


775 


681 


680 


910 


1138 


1502 




error 


15 


23 


27 


23 


19 


18 


18 


26 


30 


116 


Uf, 


ISAJET 


520 


807 


1529 


788 


714 


622 


626 


830 


1033 


1242 




SOFTSUSY 


569 


884 


1592 


866 


774 


681 


679 


914 


1142 


1297 




SPHENO 


548 


847 


1552 


828 


746 


655 


659 


880 


1080 


1266 




SUSPECT 


550 


852 


1585 


832 


754 


656 


662 


880 


1092 


1492 




error 


20 


32 


29 


32 


25 


24 


22 


35 


45 


114 


tl 


ISAJET 


379 


633 


947 


621 


523 


236 


476 


774 


951 


998 




SOFTSUSY 


398 


658 


974 


645 


544 


232 


497 


813 


987 


951 




SPHENO 


398 


658 


964 


646 


545 


248 


497 


813 


982 


986 




SUSPECT 


410 


676 


1004 


663 


560 


243 


513 


831 


1015 


1140 




error 


13 


18 


24 


17 


16 


7 


15 


24 


26 


83 



Table 5: Gluino and squark masses in GeV for the SPS benchmark points, and their 
theoretical uncertainties. The theoretical uncertainty is displayed in bold type face and is 
calculated as described in the text. 



traUno and slepton masses with accuracies at the per-cent or even per-mill level [11 
T2| , p!3| , 0, |15|. We thus take the differences in these masses accessible at -^i = 
500 GeV as a measure of the theoretical uncertainty. Tables |^ and |^ compare masses 
obtained by the four programs for the various SPS points. Only those masses kine- 
matically accessible at a 500 GeV LC are shown. The error is again defined as the 
variance of the masses as in the previous section. 

It turns out that the uncertainty in the LSP mass in mSUGRA and mGMSB 
scenarios is typically a few hundred MeV, depending on the parameter point. An 
exception is the focus point scenario (SPS2) where Sm^o = 1.4 GeV. For X2 
X^, we find uncertainties of about 1-3 GeV in mSUGRA and mGMSB scenarios. 
(Here note that with earlier program versions, especially with ISASUSY 7.51 -7.63, 
we had discrepancies of 50% and more for focus point scenarios.) For the sleptons 
we find typical uncertainties of 1-2 GeV. We note that for SPS5, the lighter stop 
would be accessible, with m^^ = 235.6, 232.3, 248.4 and 242.6 GeV for ISAJET 7.64, 
SOFTSUSY 1.71, SPHENO 2.0 and SUSPECT 2.101, respectively, corresponding to 
an error of 7 GeV, c.f. table ^ In the mAMSB scenario, SPS 9, we have much 
larger uncertainties of ~8 GeV for m^o and This is due to the fact that 

SOFTSUSY 1.71 has 2-loop GUT-scale boundary conditions for mAMSB, while the 
other programs have only 1-loop boundary conditions. If we enforce 1-loop boundary 
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conditions for all four programs the error decreases to 3 GeV. We therefore expect 
that if all programs were to use 2-loop boundary conditions the error would be better 
than 3 GeV. 

By comparing table ^ with tables we see that the present theoretical un- 
certainty in mass predictions is significantly smaller for weakly interacting sparticles 
than for strongly interacting ones. This is expected, since the former have smaller 
threshold corrections. However, the errors in tables are still larger by up to an 
order of magnitude than the expected experimental accuracies at an e~'"e~ Linear 
Collider. Moreover, the differences in the masses produce cross sections and branch- 
ing ratios that differ by a few per-cent. While this is not a problem for determining 
SUSY-breaking parameters at the low scale, it will be relevant when relating them to 
GUT-scale parameters in order to test their unification and to determine the sources 
of SUSY breaking. 



5. Conclusions 

If sparticles are detected at future colliders, measurements of their properties (and 
confirmation that they are in fact sparticles) will take place. The question following 
from this line of investigation will be: what can we learn about SUSY breaking from 
these measurements? Is there a unification of certain SUSY breaking parameters, 
and is it possible to distinguish various SUSY breaking models? 

Vital ingredients for answers to these questions will be the amount and precision 
of empirical measurements made at the particular colliders |Tl|, 0, 11^, [15], IT6|, [IS 



Another vital ingredient will be the precision with which we can relate the exper- 
imentally observed quantities at the TeV scale with the fundamental physics at 
the high-energy scale. Accurate extrapolations require multi-loop results for the 
RGEs and the related threshold corrections. In this paper, we have addressed the 
question: what is the current theoretical uncertainty associated with determining 
fundamental-scale SUSY parameters? To this end, we compared four public state- 
of-the-art MSSM spectrum calculations: ISAJET7.64, S0FTSUSY1.71, SPHENO2.0 
and SUSPECT 2.101, taking the spread of their results as a measure of the to-date 
uncertainty. Although this does not correspond to the usual notion of 'theoretical 
uncertainty', it is pragmatic in the sense that (at least) one of the available calcula- 
tional tools will be used to perform fits if and when the relevant data arrives. The 
uncertainty was shown to be largest in certain tricky corners of parameter space: 
the focus point region and high tan/?. However, even in these regions, comparison 



with previous versions of the codes |^2|] shows that the theoretical uncertainty has 
significantly improved. Sparticle masses in these regions are particularly sensitive to 
the values of the Yukawa couplings (especially the top Yukawa for the focus point, 
and the bottom Yukawa for the high tan/3 regime). Slightly different treatments of 
top and bottom masses can lead to large differences in mass predictions. It is there- 
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point 


code 


A 1 


A 2 


M ± 

Al 














la 


ISAJET 


96.6 


176.8 


176.4 


185.8 


184.9 


142.9 


201.9 


133.3 


206.0 




SOFTSUSY 


96.4 


178.2 


177.6 


188.7 


187.9 


144.9 


204.3 


136.3 


207.8 




SPHENO 


97.6 


182.9 


181.3 


190.5 


189.6 


143.9 


206.6 


134.6 


210.3 




SUSPECT 


97.4 


179.8 


179.1 


188.5 


187.5 


144.9 


204.4 


135.7 


208.1 




error 


0.3 


1.5 


1.2 


1.1 


1.1 


0.6 


1.1 


0.8 


1.0 


lb 


ISAJET 


161.2 


299.2 












195.2 






SOFTSUSY 


161.3 


303.5 












204.6 






SPHENO 


163.7 


310.2 












195.1 






SUSPECT 


162.6 


305.5 


_ 


_ 


_ 


_ 


_ 


199.7 


_ 




error 


0.7 


2.6 


_ 














2.6 





2 


ISAJET 


120.2 


221.7 


221.5 
















SOFTSUSY 


118.6 


232.0 


231.8 
















SPHENO 


123.8 


233.6 


232.1 
















SUSPECT 


123.1 


233.0 


232.6 


_ 


_ 


_ 


_ 


_ 


_ 




error 


1.4 


3.2 


3.1 














„ 


3 


ISAJET 


160.5 


297.0 








178.3 


287.0 


170.6 


289.2 




SOFTSUSY 


160.8 


300.4 








182.0 


290.4 


175.3 


292.3 




SPHENO 


162.7 


307.2 








180.2 


290.3 


172.4 


292.4 




SUSPECT 


161.7 


302.4 











182.1 


290.7 


174.8 


292.6 




error 


0.6 


2.5 











1.0 


1.7 


1.3 


1.6 


4 


ISAJET 


119.3 


219.6 


219.6 
















SOFTSUSY 


118.6 


223.6 


223.4 
















SPHENO 


121.4 


228.9 


227.4 
















SUSPECT 


121.1 


225.9 


225.7 


_ 





_ 


_ 










error 


0.8 


2.3 


1.9 














5 


ISAJET 


119.8 


225.9 


225.9 


244.3 


242.3 


191.5 


256.3 


180.9 


257.7 




SOFTSUSY 


118.7 


229.4 


229.2 


248.2 


246.1 


193.5 


259.4 


183.7 


260.7 




SPHENO 


121.1 


230.0 


229.8 


248.0 


245.9 


192.0 


259.4 


181.7 


260.9 




SUSPECT 


120.7 


230.6 


230.4 


247.8 


245.5 


193.5 


259.5 


182.8 


260.6 




error 


0.6 


1.2 


1.2 


1.1 


1.0 


0.6 


0.9 


0.7 


0.9 



Table 6: Differences in predicted masses in GeV for the SPS points la to 5 (mSUGRA 
points). Only sparticles with masses that kinematically can be produced at a 500 GeV 
Linear CoUider are displayed. The theoretical uncertainty is displayed in bold type face 
and is calculated as described in the text. 



fore critical that the accuracy in deriving the running top and bottom DR Yukawa 
couplings is at a maximum in any calculation. 

We used previous LHC estimates of expected empirical errors at two benchmark 
points to perform separate fits to mSUGRA with the four different codes. The 
parameters resulting from the four calculations show a difference comparable to the 
statistical error upon them, showing that theoretical uncertainties must be taken 
into account. We then went on to quantify the theory uncertainty associated with 
squark and gluino masses for the SPS benchmark points. The theory uncertainty 
was also quantified for parts of the MSSM spectrum that would be kinematically 
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point 


code 


A. 1 


A.2 


M ± 














6 


ISAJET 


190.0 


218.1 


215.7 






236.8 




227.8 






SOFTSUSY 


190.2 


220.6 


217.9 






240.9 




232.8 






SPHENO 


191.2 


225.3 


222.3 






237.6 




229.1 






SUSPECT 


190.2 


222.2 


219.6 


_ 





240.9 





232.2 







error 


0.3 


1.7 


1.6 






1.3 




1.4 




7 


ISAJET 


162.4 


268.0 




248.7 


248.3 


127.3 


261.1 


119.9 


263.3 




SOFTSUSY 


163.6 


263.5 




247.2 


246.9 


126.4 


259.3 


120.5 


261.2 




SPHENO 


163.4 


271.1 




251.6 


251.3 


131.0 


265.3 


123.8 


267.3 




SUSPECT 


163.6 


262.2 




246.9 


246.6 


127.8 


259.4 


121.6 


261.4 




error 


0.3 


2.4 




1.2 


1.2 


1.2 


1.6 


1.0 


1.6 


8 


ISAJET 


137.4 


254.6 








175.7 




168.9 






SOFTSUSY 


138.4 


261.2 








175.4 




169.9 






SPHENO 


139.2 


266.1 











180.3 


_ 


173.5 







SUSPECT 


140.0 


263.5 


_ 


_ 


_ 


177.6 


_ 


171.8 


_ 




error 


0.6 


2.8 








1.3 




1.2 




9 


ISAJET 


174.8 




175.0 
















SOFTSUSY 


196.7 




196.7 
















SPHENO 


168.0 




168.4 
















SUSPECT 


167.3 




167.3 
















error 


7.9 




7.9 















Table 7: Differences in predicted masses in GeV for the SPS points 6 to 9 (6: non-minimal 
SUGRA, 7+8: mGMSB, 9: mAMSB). Only sparticles with masses that kinematically can 
be produced at a 500 GeV Linear Collider are displayed. The theoretical uncertainty is 
displayed in bold type face and is calculated as described in the text. 



accessible to a 500 GeV Linear Collider for the SPS points. The tables of spectra 
listed in the present article may be used for further comparisons with other (possibly 
new or private) computations. The theoretical uncertainties for LHC and Linear 
Collider observables were shown to significantly depend upon the point in parameter 
space being considered. It will not be very precise, therefore, to use the theory 
uncertainties calculated here when analysing actual data, since the SUSY breaking 
parameter point is a priori unknown. 

If sparticles are measured, the correct modus operandi is clear when analysing 
data: any fits to a fundamental SUSY breaking model should be performed with 
several state-of-the-art spectrum calculations in order to deduce the level of theo- 
retical uncertainty, which must then be included in final quoted errors. The theory 
errors must clearly be included when discussing the accuracy required to distinguish 
between SUSY models on sparticle observables, such as in reL . 

We note that each of the codes used assumed a pure MSSM desert up until the 
GUT scale ~ 10^^ GeV. We know already ||51, 52] that this assumption is called 
into question because the gauge couplings themselves require ~ C(l%) corrections 
at the high scale in order to make them unify properly. It is thought that this 
level of correction is perfectly reasonable in, for example, SUSY GUTs since model- 
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dependent GUT-scale threshold corrections |53, 54] are expected from (for example) 
heavy coloured triplets. If the combined theoretical and empirical accuracy is signif- 
icantly better than 1%, then observables at colliders could be used to measure these 
threshold corrections. 

Fortunately, theoretical errors in sparticle mass predictions are certainly not 
static. There has been much progress reducing them recently (especially in the Higgs 
and electroweak symmetry breaking sectors, see for example [^, ^ |55|, 0). We 
expect this trend to continue, which, as our present results indicate, is desirable if 
we are to disentangle SUSY breaking from experimental observables. 
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